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ABSTRACT 

Multiloop response properties of controllers are. In general, very 
difficult to obtain because an independent forcing function is needed for 
each describing function to be measured, and interpolation procedures may 
be required to obtain intermediate describing functions at common frequen- 
cies. Even then, a certain amount of untangling is required before the 
final results are obtained. When the loops which ire closed and the nature 
of *.he describing function forms adopted in each loop are known or hypothe- 
sized, matters can v .* made . -ch sir-ler. Then, the quantitative values of 
the individual describing functions readily be identified using appro- 
priate closed-loop describing function measures and decomposition procedures. 
Two examples are provided for the measurement of driver /vehicle multiloop 
response properties using a single disturbance input. The validity of the 
procedure is bawl on current multiloop operator adjustment rules and is 
made plausible by comparison with experimental data. 


BfTHQDUCTIOH 


The mo*t common nun/machlne system in use today is a driver and an 
automobile, yet remarkably little is known empirically about the details 
of driver dynamic responses and how these Interact with the vehicle dynamic 
characteristics. For instance, only one study of lateral control has been 
reported where driver describing functions were measured in an actual car 
'Ref. ] ), and less than half a dozen studies have even been concerned with 
measurements in a simulator situation. A primary reason for this lack of 
attention resides in measurement difficulties due to the multilOQp nature of 
most driver/vehicle system controx situations of interest . The difficulties 
item from the need for one independent input for eaeft describing function to 
r-' measured, to the many couplings te tween the vehicle's state variables, 
to jignal-to-nolse problems, and last but not least, to data manipulation 
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difficulties which include fairing and Interpolation with already low confi- 
dence, dubious data and snail differences between large uncertain lumbers 
(Ref. 2), In spite of these difficulties, enough measurements now exist to 
support sons key rules for multiloop human operations (Refs. 2-U). When 
these rules can be applied to specific situations, the fundamental diffi- 
culties can be bypassed and much simpler measurements can be made which 
nonetheless reflect the slgnf leant multiloop human control characteristics 
(Refs. 5-6). In this paper these techniques art applied to the problem of 
driver steering control. 

VSOGLg A® DRIVER RMTOBM DYMMIC8 

The object of control in the driver /vehicle system is the automobile. 

Its dynamic characteristics for essentially constant speed steering maneuvers 
can adequately be described by equations of motion involving side velocity, v, 
heading rate, r, and roll angle, <p. The roll mode is least Important in «*wige i 
maneuvers and many of the roll effects, such as roll steer and camber thrust, 
can be partially accounted for in the lateral velocity and beading degrees of 
freedom. Accordingly, the constant velocity lateral-directional character- 
istics can often adequately be approximated by the two-degree-of- freedom sett 



Approximate values for the stability derivatives in terms of vehicle para- 
meters and the key vehicle transfer functions for control inputs are eumarlted 
in the appendix. More complete automobile descriptions in six degrees of free- 
dom are provided in Ref. 7. 

Turn now to the driver. Driving tasks are often multiloop in nature. 

That is, the driver responds to more than one vehicle motion quantity. The 
key to multiloop driver models is the fundamental eoneept of manual vehicular 
control analysis: that tbs operator constructs feedback loops about the effec- 
tive controlled element. The feedback quantities available to him for possible 
use consist of those: 
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• Directly sensed within the general visual field. 

• Observable via visual displays. 

• Directly sensed using modalities other than vision. 

^iuantities which can be perceived from the fixated visual field will show 
no scanning penalties, whereas those which require Instrument scan or modi- 
fi cation of the fixation point will introduce decrements in some driver 
dynamic features. For lateral steering control the feedbacks are derived 
from the general visual field and possibly supplemented by non* visual 
modalities. The feedback quantities actually selected by the driver will 
be those needed to aatj.sfy the guidance and control needs and certain 
driver-centered requirements. The guidance and control needs are sltuatlon- 
3oecific. In the steering case, they involve an outer loop for lateral 
position control and an appropriate inner loop to provide damping for the 
low-frequency mode formed by virtue of closing the lateral position loop. 
These requirements are, of course, essentially Independent of whether the 
controller is animate or inanimate. 

The driver-centered requirements are central to the manual control, as 
opposed to the general control, problem. The human propensities and behavior 
associated with these requirements must be discovered by experiment. From 
the data available (Refs. 2 - 4 ), a series of adjustment rules similar to those 
for single-loop manual control systems can be stated. These include: 

i . The feedback loops preferred are those jfhlch: 

a. Can be closed ’ th minimum operator equalization, 

b. Require minimum scanning. 

c. Permit wide latitude in variation In the operator's 
characteristics . 

T. Where distinct Inner- and outer-loop closures can be 
defined by ordering the bandwidths (i.e., the higher 
the bandwidth, the more inner the loop), a series 
multiloop structure applies. 

v . The adjustment of the variable gains in each of the 
loops is, in general, similar to that which would be 
used by a skilled automatic control designer who has 
available the snme feedback entitles. To a first 
approximation: 

a. The crossover model Is directly applicable to many 
inner loop closures. 
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b. The crossover model also holds for the outer loop 
with the proviso that the effective cont rolled- 
element transfer function includes the effects of 
all the inner-loop closures. 

4, When scanning is not present, the remnant Is primarily 
associated with the inner loop and is essentially the 
same as that for a single -loop system squi valent to 
the inner loop alone. 

Other rules apply to situations where scanning is present } these are not 
pertinent hers. 

rmm/maasM smm bt h uo t u bs 

The combination of driver and vehicle into an appropriate control system 
for lateral position can conceivably be accomplished using a wide variety of 
feedback quantities. A number of these have been Investigated theoretically 
in Refs. 6 and 9 . One of the most likely possibilities when guidance and 
control requirements, driver -esnttred requirements, availability of cuss ir 
the visual field, and interpretation of such experimental evidence as driver 
eye movements are considered Is an outer loop in lateral position, y, and an 
inner loop Involving either path angle or heading angle. Additionally, path 
rate and heading rate ars pertinent for higher- frequency control action. 

When all evidence is taken into account, * very likely structure for many 
driving tasks is that shown In Fig. 1 . In squat ion form, the quasUlneer 
model for the driver is: 

- Yy(ja>)Yf(Jo)y t (Ja>) - Ju>)*( Ja>) ♦ Tf(Ju>)n( Jc») (2) 

**Jhs describing fi action, Y y Y^, characterizes the driver's operations on 
functions of the position error, whereas Yf alone represents ths driver's 
operations on functions of heading angle. The portion of the driver's 
output which is not linearly correlated with any disturbance or coanand 
inputs is modeled by ths remnant, n. 

It is important to note that the above presumed structure in no way 
ImpHas direct perception of either lateral position or heading angle as 
such, but only that th* driver responds in part to some function of these 
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>ck Diagram of Driver- Vehicle 
Lateral Steering Control 


variable* a* sensed from the available visual field and non-visual madallties. 
In other word*, although Bq. 2 it expressed explicitly in tens* of lateral 
position and heading angle, any other possible feedback* which ere linear 
function* of y end q arc aleo handled Implicitly. 

When y end f actually are the direct function* used by the driver in 
th* experimental situation of Intcreet, then in the light of the oultiloop 
principle* listed abeve heading should be an inner and lateral position an 
outer loop, and the low- end mid- frequency characteristic* of Yy should be 
a simple gain and those of Yq a gain plus tins daisy, dines sens vehicle 
dynanic* present c ha lle n ge* which ere readily overcome with handle rata, 
the higher-frequency portion of Yq should include a lend. This can be con- 
sidered as either high-frequency lead generation on ♦ (with little or no 
degradation in effective time delay) or direct sensing of yawing velocity. 

The latter lc pertinent for supra- three hold values of yawing velocity, r, a* 
sensed by the vestibular system , Thus, the driver dynamics which we will 
ultimately use arc given tyi 


*y " *y 

Yq - ♦ O 


(3) 


Mof* extensive high-frequency characterisation can bo used with appropriate 
changes in r. Per example, when distinctions ere to be drawn between different 
power steering systems, or when the effect* of alcohol on the neuromuscular 
eyetem arc to be considered, a more extensive analytical description of the 
neuromuscular system cam be added to Y, with t adjusted accordingly. 

If interpretation of data obtained using these specific forme are in 
reasonable accord with analytical implications, then this would be strong 
Indirect evidence that the driver used q and y feedbacks directly In aceeop- 
llehlng the particular driving task considered. If, on the other hand, Y y 
or Yq as deduced from the task-specific data are much more complex dynamic 
forma, then alternate feedback* should be considered. 
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In the study of driver/ vehicle interaction in steering tasks, a number 
of system inputs require consideration. The most obvious ere the lateral 
positional command, y c , and external disturbances, T), because they are 
representative of actual roadway conditions. As shown in Fig. 1, y c Is 
actually present in the driver component of the closed-loop system, i.e., 
the lateral position command actually acting on the system is to some extent 
driver induced. The general character of the command is, however, determined 
bv the roadway, lane pattern, obstacles to be avoided, etc. It can be as simple 
as a prescribed pathway with narrow tolerances and as complex as a tortuous 
way through freeway traffic. The external disturbances can arise from the 
atmosphere, such as yawing velocity, r g , or side velocity, v g , gusts, from 
roadway- induced disturbances, or from specially-contrived force and/or moment 
generators attached to th. vehicle (e.g., a rocket). In addition to these 
inputs into the driver/vehicle system, various test disturbance are of interest 
for special measurement applications. The three most common are indicated in 
Fig. i . The first, & d , is a front wheel steer angle applied in series with 
the driver's steering wheel input. The 6 d input Is readily applied In actual 
or simulated cars by the addition of an extensible link or other differential 
device into the steering system. The heading and path disturbances, and y^, 
are primarily of values in simulator applications, whe"* these signals are 
readily added within the equations of motion. These test disturbances enter 
the equations of motion at the locations shown In Fig. 2. They can, for 
example, be readily applied to the servo drive of a television camera used 
to generate the visual scene from a model landscape, as in Refs. 10 and 11 . 

The equations of motion for lateral position and heading with these forcing 
functions and disturbances and using Yy and for generality are given by: 

(i *r,w4) |Y1 «V» < W*°S V«°? a? 

■ 

W»°6 (’♦wb * 0 1 4 W«°J 
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Figure 2. Test Disturbance Entry Locations 
In Vehicle Dynamics 


Ths front wheel tteer angle is provided by the auxiliary equation: 

ft - S 4 ♦ y) + YyO^n - *) (?) 

The closed-loop system responses of heading, lateral position, and front 
wheel steer angle are given as functions of the forcing function and dis- 
turbances in Table 1 , The multUoop character of this system is indicated 
by the sum of Y^GgGjjj and YyY^OgO^ in the denominator and the presence of 
the coupling numerator terms, and Q*g, shown in the * and y numerator 
expressions (Ref. 12). The two prises on the denominator indicate that two 
loops are closed. 

Closed-loop describing Auction measurements can be made using any of 
the forcing Auctions or disturbances, or combinations of these which are 
independent. Because the remnant acts as a continuous power spectral den- 
sity, it will tend to contaminate any such measurements as a noise* However, 
If input power is large * compared with the remnant when both are reflected 
to a point of interest in the system, this contamination will be negligible. 

In any event, it can generally i assessed directly by consideration of the 
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TABLE 1 

CLOSED-LOOP SYSTEM DYNAMICS 


Denominators 


D" - 1 ♦»,(»,«£♦ YyOg) - 1 + + 


y Numerator; 


i 

& 

© 

*- 

i 


D-y - V*G,agy c + (l + Y + 0,oJ)y d + Y^ofr + 0&, ♦ (0^ + Y^O^Jr, ♦ (£ - Y^fl, ^ 


* Numerator: 


D "+ * Vft 0 -* - V*°«°» y d + We" + + < a J + V> 0 8 °$ )T1 + (’ + V*°8 $)*d 


6 Nun.^rator: 


D"6 . Y y Y t G,y c - Y^O^ 4 Y^J.n + & d - V.^ +0 P’> “ Va(' *^)»d 



relative power at the frequencies Involved. One of the tricks In mult Hoop 
measurement is to select a family of describing functions which indicate 
high Sigiial-to-noise ratios over particula* (different) frequency ranges. 

The fundamental rscaj-^rement problem is to determine y y» Y^, and remnant 
po v<r spectral densities, In principle, to obtain Yy and Yy separately, 

two independent inputs will be required. As noted already, the early work on 
suit Hoop human operator measurements had to use several independent inputs 
because the general nature of quantities such as Yy and Y^ were to be dis- 
covered. How that we have some appreciation for at least the likely form of 
these quantities, as given in Eq. 3, a simpler procedure can be used with only 
one input. The Inputs to be considered she id satisfy two fundamental criteria: 

• The cicsed-loop describing functions should have a large 
frequency range over which signal to noj.Ee ic high, i.e., 
remnant power should be small compered to forcing function 
or disturbance power at common points with the system. 

• The closed-loop describing functions should be differ- 
entially sensitive to Y y and Y^ in different frequency 
bands so that the properties of Yy ±nd Y^ can be readily 
untangled. 

For these reasons, the principal disturbances and forcing functions that have 
been used are y c , 5 d , * d , and y d . fee first two are applicable both actual 
vehicles on the road and simulator operations, whereas the latter two are 
suitable only in a simulator. Measurements with a disturbance, while having 
grtiat value in peculiar situations, tend to be extraordinarily muodled because 
of the coupling features introduced via the coupling numerators. It is better 
to stay with inputs acting directly within one of the loops rather than a 
disturbance which is diffused into both loops viu complex vehicular couplings. 

Two particular examples will be illustrated in more detail below. These 
arc describing functions measured with a steering or with a heading disturbance 
as the test input. 

TYPI CAL SR0CHXBIHQ OTCTJ0JI8 KISH 0KKBR33 OI0TURBAKXS 

The front wheel steer angle response with a steering disturbance as the 
ays* -n forcing characteristic and remnant is given by specialisation of 
Table 1, i.e.: 
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Here the quantity Yj is seen to be obtained by taking the closed-loop 
response measurement, 6/d d , inverting it, subtracting dividing out the 
known transfer function of the steering linkage and vehicle, OgQ$. YJ la 
particularly simple In fora in that the driver’s b ea din g describing function, 
Yf, is a multiplicative factor, while Yy enters into the bracketed expression 
in a relatively simple way. 

At this point, let us examine the properties of the bracketed expressions 
In Eq. 8, In terms of the lateral equations of ths automobile, t’lis Is given 
by: 
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<7,u, |s| (Log Scot#) 

a ) Conventional Root Locus b) Bods Root Locus 


Figure 3. 


System Survey of Hoots of ) ♦ Ky{o£/oJ[) 



Using the relationships developed in the appendix for the two-degrees of 
freedom equations this 1 ecomes: 



The roots of Eq. 9 can be readily evaluated if it is viewed as the result 
of a loop-closing operation in which the transfer function. Yy j f, is the 
open-loop transfer function. At this point, vc bring to bear the multiloop 
operator adjustment rule*; previously described and subat*tute for the general 
outer- loop driver describing function, Y y , the much more specific form, Ky. 

A system survey which shows tbe roots of 1 + M \ is given in Fig. 3. The 
Bode root locus there indicates the closed-bxip roots as the driver gain, K^, 
is varied. Tht conventional root locus, which shows the same effects, is 
remarkable in that it is an example from practice of an academic anomaly, 
i.e., the damping of the quadratic zeros is approximately the same as the 
heading transfer function numerator inverse time constant, l/T r . This 
approximation is exact for the two-degree- cf- freedom lateral-directional 
equations when the z axis radius of gyration, k z , Is equal to the geometric 
mean between the car dimensions, a and b. 

The path loop lb the outer loop, so typically the gain, Ky, is relatively 
low. Then, a nominal zero dB line on the Bode root locus would be as shown 
in Fig. 3 . The a -Bode (locus of all real roots) is essentially on the low- 
frequency asymptote at this point; consequently, one of the roots will be 
given almost exactly by: 


*1 s -Uo*y 


00 


Because of the peculiar approximate relationship between 2(&cu)y and 1 /T r , 
the sun of the closed-loop roots is constant in spite of the variation in 
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gain (this is not, of course, true for a general open-loop system with a 
quadratic numerator and denominator). Under these special circumstances, 
the second root will be given by: 


Cg * — { l/Py " TT -»Ky) 

- 


With these literal approximate factors, the value of 1 + Ky(oj£/ojt|) becomes: 




(13) 


When the forer for Y^ provided by the multJlocp adjustment rules is used in 
conjunction with Eq. 13, the total describing function, Yp, is seen to have 
the form: 


+ 1 


+ i)(t^ + 


An asymptotic Bode plot of this is shown in Fig. A. Here, it will be noted 
that the low-frequency asymptote defines the outer-loop gain uniquely, while 
the mid- frequency horizontal asymptote Is Just tbe gain, K*. The breakpoint 
between the very low-frequency and mid- frequency asymptotes is U^. Tbe 
high- frequency breakpoint gives the indication of any inner-loop lead. Actual 
data are presented on Fig. 5* They Bhow the same characteristic concave- 
upward basket-llke shape of the sketch of Fig. 4. Quick approximations to 
the values of Ky, K^, and 1 /? L are readily obtained by fairing in the asymptotic 
plot. By keeping in mind the Shape and functional relations involved in the 
asymptotes, comparative describing function data and parameters can be directly 
assessed by eyeball. 

Precise values for the driver characteristics are most easily determined 
with an optimal curve-fitting routine. The routine adjusts the Yp parameters 
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including the effective time delay, t, to simultaaeoualy fit amplitude and 
phase* Ad exwnple curve fit Is also mmu on Fig. 9 (although the T r , Tr 
dipole in Bq. lU was not included for simplicity). 

As a practical matter, the measurements are quite simple to obtain with 
sums of sinusoids and on-line describing function measurements. These are 
then inexpensively reduced and curve-fitted on a routine basis. Thus, the 
previously highly- complex multiloop reduction procedures are reduced to an 
elementary and inexpensive set of operations. The shape of the data further 
confirms the multiloop describing rules presented earlier, in that the data 
behave according to the predictions. 

Yj data can be interpreted in alternate ways which are helpful When 
considering actual perceptual structures which the driver may he using. For 
example, excellent single-loop driver/vehicle systems result when the driver 
Is assumed to at' on the “time-advanced lateral deviation, 0 y(t+T), or the 
"aim- point-heading" error, Both are shewn in Fig. 6 and both Involve an 

aim point at a distance, B, down the road. With the first: 


y(t + T) * y + vr 
» (Te + l)y 


(15) 





In this interpretation, the low-frequency lead breakpoint in Fig. U is Just 
1 / T. Thus: 


1/T - 

UoKy I U 0 /R 

06) 

*y 1 

1/R 

(17) 


the inverse of the aim point distance down the road. With the aim point 
heading error: 

* L = + + y/R 08 ) 

If the composite feedback quantity, Y^(* + Kyy), in Eq. 2 is identified with 
? L , then Ky i 3 simply l/R, Ju3t as above. With either of these variables as 
the basis for control, there is an effectively single-loop closure, with the 
kinematic relation l/R corresponding to Ky, and ] or R/U 0 to the low- 
frequency lead time constant. 

MULTILOOP INSCRIBING FUNCTIONS 
WITH HEADING DISTURBANCE 

The closed-loop heading-to-heading-disturbance transfer function as 
obtained from Table 1 is: 


(19) 


When this . considered in a fashion similar to that above, i.e., when 
is taken as an input and * as the error of an equivalent single-loop feed- 
back system, then the effective inner-loop describing function, (t d />Sr) - 1, 
becomes: 
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Here, the effective vehicle transfer function, [G^]', is the transfer function 
for the vehicle as seen by the heading loop as an outer loop, i.e., it reflects 
the lateral position loop closure. If, now, we insert the appropriate forms 
for vehicle anJ driver multiloop characteristics, Eq. 2C becomes: 


w»(i fe+1( Tt> + 1 ) e ~* OTt 


An asymptotic Bode plot ox' this function is shown in Fig. 7. Here again the 
path loop dominates the low- frequency characteristics, while the heading 
inner loop Is more ilgnif leant at higher frequencies. In Fig. 7 note that 
the pole at very low- frequencies is non -minimum phase; hence, the phase curve 
starts at -270 deg. This pole arises from the terms in brackets in Eq. 21 . 
This function can be used directly to fit data, as with the steering input 
considered previously. Typical data from Ref. 11, as well as same unpublished 
data from that series of experiments, is shown in Fig. 8. These data indicate 
a value of tyCy of about 0.37 rad/sec whereby the low-frequency lead time con- 
stant is about 2.6 sec and, for a speed of 88 ft/sec (60 arph), the "aim point 
distance" would be about 2^0 ft. 
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upper middle to hi^h frequencies, the effective vehicle transfer 
fur. *lon, [C? ' , is hardly affected by the path closure. Thus, at these 

froq * ^ncies: 



The driver's heading response describing function, can then be isolated 
fr_*a tne equivalent open-loop system high-frequency properties by removing 

d/ 

th*-» '.nown vehicle and steering subsystem dynamic characteristics. G a 0£, 

Typical data, from Ref. % J, are snown in Fig. 9. These demonstrate* clearly 
the difference between Subject B and Subject F ir that only the matter uses 
nigh- frequency lead generation. Again, th^* fat that data, tueh as shown 
in Figs. 1:5 and 9, tend to follow the trends predicted using the multiloop 
adjustment rules provides further verification. 

Finally, an interesting comparison of the features of describing functions 
oc'tPined from heading vs. steering disturbances can be seen in the low-frequency 
regions of Figs, t and 8. The Yp steering disturbance case has effects 

changing the low- frequency amplitude while the phase is insensitive (etarting 
at —90 deg for any value of Ky). Just ttio opposite occurs Cor the heading 
disturbance case in Fig. 8 where the amplitude is relatively insensitive, 
while the phase curve moves directly with Ky. 

CONCLUSION? 

The point of view taken in thi3 paper permits the simple and routine 
estimation of the mere significant driver characteristics in steering control 
ta Ihc fundamental difficulties of multiloop measurements are avoided by 

tailoring the driver loop structure per the currently understood mult iloop 
human operator adjustment rules. The simplifications so obtained permit 
Tulttl.iop ^ar/driver situations to be analysed as easily as systems which 
ar« -a sent tally single loop and single Input. 
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NOMENCLATURE 

Distance of mass center aft cf the front axle 
Distance of mass center forward of the rear axle 
Steering linkage transfer function 

Vehicle transfer function between steer angle and lateral 
velocity ' 

Vehicle transfer function beeween steer angle and heading 
angle 

Vehicle yaw moment of inertia 
Vehicle radius of gyration, Jl zz /n 
Driver gain for lateral position control 
Driver gai.j for heading control 
Vehicle mass 
Driver remnant 

Ya-./ing velocity (heading angle rate) 

Lapluca transform complex variable 
Nominal forward velocity 
Side velocity 

Lateral position of mass center (relative to an inertial 
frame) 

Lateral position cowaand to driver 
Lateral position disturbance 

Driver describing function for lateral position control 
Side force due to front tl~e slip angle 
Side force due to rear tire slip angle 
Driver describing function for heading angle control 
Path angle 

Front wheel steer angle 
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Steer angle disturbance 
Driver steering wheel >vngie 
General external disturbance 
Effective time delay 
Heading angle 

Heading angle disturbance input 
rrequency 
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The two-degree-of-freedom latf ral -direct ional equations for a car are 
given in Eq. 1 . The stability derivatives are defined in terms of vehicle 
and tire design parasieters by the expressions below. 
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ABSTRACT 


AK EXPERIMENTAL STUDY OF THE MOTORCYC'E 
ROLL STABILIZATION TASK 

David J. Eaton 

The stabilization task of tha motorcycle rider ha* been studied 
experimentally, using data from actual road tests, In each experi 
ment, the test subject, whose upper body was restrained by a rigid 
brace, controlled the motorcycle roll angle by applying a steering 
torque to the handlebars. Roll angle, roll rate, steering torque, 
and steering angle were recorded for three test subjects. Experi- 
ments ere restricted to the s tcady -s tate , straight path, constant 
forward speed (usually 30 mph) situation. 

The strongest source of excitation to the man-motorcycle 
system was found to be the rider's remnant. To minimize bias errors 
in identifying a linear transfer function representing the rider, 

8 technique developed by hingrove and Edwsrds of the NA! < 
Research Center was employed. The result! are consistent with a 
theoretical analysis of the man-motorcycl a system, using classical 
automatic control theory. 



